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ABSTRACT 

A multi-transition 3 mm molecular line single-pointing and mapping survey was carried 
out towards 29 massive star-forming cores in order to search for the signature of inward 
motions. Up to seven different transitions, optically thick lines HCO"'"(1-0), CS(2-1), 
HNC(l-O), HCN(l-O), 12CO(1-0) and optically thin lines Ci«O(l-0), i3cO(l-0) were 
observed towards each source. The normalized velocity differences (i5VcSj i5Vhco+) 
between the peak velocities of optically thick lines and optically thin line C^^O(l-O) 
for each source were derived. Prominent inward motions are probably present in either 
HCO"'"(1-0) or CS(2-1) or HNC(l-O) observations in most sources. Our observations 
show that there is a significant difference in the incidence of blue shifted line asymmet- 
ric line profiles between CS(2-1) and HCO+(1-0). The HCO+(1-0) shows the highest 
occurrence of obvious asymmetric feature, perhaps owing to different optical depth 
between CS(2-1) and HCO'^(l-O). HCO+(1-0) appears to be the best inward motion 
tracer. The mapping observations of multiple line transitions enable us to identify six 
strong infaU candidates G123.07-6.31, W75(OH), S235N, CEP-A, W3(0H), NGC7538. 
The infall signature is extended up to a linear scale > 0.2 pc. 

Key words: ISM: molecules - ISM: kinematics and dynamics - stars: formation - 
radio lines: ISM . 



1 INTRODUCTION 

, Massive stars, once formed, dominate the luminous, kine- 
' matic, and chemical output of stars and their feedback play 
a dominant role in the evolution of molecular clouds and 
any subsequent star formation therein (Krumholz & Bon- 
nell 2007). Despite their importance the formation of mas- 
sive stars is currently an unresolved problem. The difficul- 
ties lie in their short evolution time scale and complex star- 
forming environment in star cluster, typically embedded in 
the dense cores. Two main theories for high-mass star for- 
mation are mergers of several low-mass stars (Bonnell et al. 
1998), which requires a high stellar density; and accretion of 
materials from the circumstellar disks, which requires high 
accretion rate to resist radiation pressure. The latter is sim- 
ilar to that of low-mass star formation (Shu et al. 1987), 
which accompanied by outflow during the process of grav- 
itational collapse. It appears that more observations favor 
the latter scenario for massive star formation (e.g. Zhang et 
al. 2001; Beuther et al. 2002; Patel et al. 2003; Moscadelh 
et al. 2005; Jiang et al. 2005; Fuller et al. 2005; Behran et 
al. 2006). Nevertheless, intermediate theories that involve 
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both accretions and mergers are now also quite commonly 
accepted (e.g. Whitney 2005; Bonneh & Bate 2006). 

While accretion often can not be observed directly, its 
presence can be inferred from the presence of large scale in- 
fall and outflow (e.g., Klaassen & Wilson 2007). Therefore 
the infall motion is one of the important elements to help 
understand the theories of massive star formation. These 
motions can be studied by investigating the profiles of opti- 
cally thick molecular lines, which show the blue asymmetric 
structure, named "blue profile", a combination of a dou- 
ble peak with a brighter blue peak or a skewed single blue 
peak in optically thick lines (Mardones et al. 1997). Before 
a definite claim of collapse signature from blue profile, we 
should rule out all other possibilities. For example, we know 
double peak could be caused by two velocity components 
in the line of sight. To rule out this possibility we must ob- 
serve optically thin line that peaks at the self-absorption dip 
of optically thick line (Gregersen et al. 2000; Wu & Evans 
2003). Surely blue profile may also be caused by rotation 
or outflow. However, infall motion is the only process that 
would produce consistently the blue profile. Outflow and ro- 
tation only produce a blue asymmetric line proflle along a 
particular line of sight to a source. Outflows strongly affect 
the line wings in the optically thick molecular line along the 
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outflow axes. Most infaU candidates both show infaU signa- 
tures and young outflows, but outflow docs not seem to be 
responsible for blue asymmetry alone (Klaassen & Wilson 
2007). Rotation might cause blue asymmetric line profile at 
one side of the rotation axis, but at the same time the red 
asymmetry could appear at the other side of the rotation 
axis. Mapping observation allows us to locate the center of 
the inflow motion and to identify dense cores that are simul- 
taneously showing evidence for inflow motions and outflows 
(Wu ct al. 2007). 

Blue profile was detected in B335, a low-mgiss star 
formation region (Zhou et al. 1993, 1994), which is well 
matched by the standard model of "insido-out" gravitational 
collapse (Shu 1977). Subsequently, more and more blue pro- 
files were found via large surveys towards both young stel- 
lar objects (Mardones et al. 1997; Gregersen et al. 2000) 
and starless cores (Lee et al. 1999, 2001, 2004; Sohn et al. 
2007) in low-mass star formation regions. If massive; stars are 
formed by accretion/collapse in a similar way as low-mass 
stars, blue profile of optically-thick lines should also be found 
towards them. Actually, these inflow signatures have been 
reported in many massive sources (Welch et al. 1988; Dickel 
& Auer 1994; Wolf-Chase & Gregersen 1997; Zhang, Ho & 
Ohashi 1998). Especially, several large surveys of inflow sig- 
natures in optical-thick lines towards massive young objects 
have been carried out, which give the important evidence 
of inflow/accretion/coUapsc (Wu & Evans 2003; FYiller et al. 
2005; Klaassen & Wilson 2007; Wu et al. 2007). In these sur- 
veys, the infall velocity is found to be between 0.1 kms~^ 
and 1 kms~^ and the mass infall rate is estimated to be 
between 2xlO"''M©/yr and 10"^M©/yr (e.g. Puller et al. 
2005). 

In this paper, we identify and study the inflow mo- 
tion using several optical-thick lines including HGO"'"(1-0), 
CS(2-1), HNC(l-O), HCN(l-O), and ^^00(1-0). Our inves- 
tigation of the profiles of optically thick molecular lines is 
aimed to understand two questions: (1) Is the infiow mo- 
tion in massive star formation cores similar to that in low- 
mass star formation cores? (2) Which optical thick molecu- 
lar line is the best tracer for the infall motions in massive 
star-forming cores? Actually, such mapping observations to- 
wards collapsing massive star-forming regions, yet are sel- 
dom, enable us to identify the most probably infall candi- 
dates. Therefore, our observational study of multiple lines in 
29 high-mass star-forming regions can help better character- 
ize the detailed properties of massive star formation in the 
dense cores, reveal more information about the collapsing 
signatures and provide some strong infall candidates. 

After details of our sample selection and the related ob- 
servations in the next section, the main results of this paper 
are presented in §3; §3.1 gives the single pointing observa- 
tion results; §3.2 presents the mapping results and analyses 
towards the infall candidates. We discuss our observational 
results in §4 and summarize in §5. 



2 SAMPLE &: OBSERVATIONS 

As a whole, the sources were selected by applying the fol- 
lowing criteria: (1) they are associated with water masers 
(24 sources) or massive young stellar objects (5 sources); 
(2) they were included in the Spitzer Space Telescope pro- 



grams before 2007; (3) source distance is less than 4kpc. All 
of the 29 sources in our sample are listed in Table 1 with 
22 sources (sources above the blank line in Table 1) selected 
from Shirley et al. 2003. To enlarge the sample, we also se- 
lected 7 sources (sources below the blank line in Table 1) 
from Spitzer programs that have images of massive star- 
forming regions. Among the 7 sources, S88 and IRAS19410 
are also listed in Shirley et al. (2003) with different names as 
S88B and G59. 78-1-0. 06 respectively, but the source center 
positions are slightly different. 

There is evidence that water masers axe a signature of 
recent or ongoing high-mass star formation (e.g. Tofani et al. 
1995). The first criterion of our sample selections can ensure 
that all sources in our sample arc associated with massive 
star-forming regions. With high sensitivity and high spa- 
tial resolution at infrared wavelengths, the Infrared Array 
Camera (IRAC; Fazio et al. 2004) and Multiband Imag- 
ing Photometer for Spitzer (MIPS; Rioke et al. 2004) see 
through dust into the birthplaces of stars and provide new 
insight to directly reveal the distribution of star formation 
sites and their environments in which stars form. IRAC data 
are available to all of the tragets and the MIPS data are also 
available for most of thcni. Considering the limited resolv- 
ing power of the telescopes and the size of the infall region, 
we made some efforts to exclude distant sources. The source 
distances in the resulting sample are in a range from 0.6 to 
4.0 kpc. The positions and distances are mainly determined 
from two literature search (see Table 1 for references) . Most 
of the sample include cores of ultra-compact Hii (UC Hii) 
regions, the sample also include cores of compact Hii (CHii) 
and Hii regions. 

Our observations were carried out using the 14 m tele- 
scope of Purple Mountain Observatory (PMO) in Delingha, 
China. The position-switch mode was used, the pointing 
accuracy is estimated to be better than 9". One SIS re- 
ceivers operating at ~110 GHz were used in the observa- 
tions along with an acoustic-optic spectrometer (AOS) back 
end with 1024 channels. The standard chopper wheel cali- 
bration was used during the observations to get the antenna 
temperature, TJ, which has been corrected for atmospheric 
absorption. The typical system temperatures Tsys is about 
250-350 K. At central source position, the accumulated inte- 
gration time towards each source is larger than 10 minutes 
and the signal-to-noise ratio is larger than 10. We carried 
out our mapping observation only towards strong emission 
(Ta > 0.5 K in HCO+(1-0) and CS(2-1)) sources with line 
asymmetry for the sake of the observational efficiency. The 
positions of CO emission peak (which were marked in Ta- 
ble 1) have the highest priority in further HCO"'"(1-0) and 
CS(2-1) observations. 

The mapping step is 30" for CS(2-1), HCO+(1-0), and 
1' (~ beam size, for our source distance, corresponding to 
linear resolution between 0.16 pc and 1.07 pc) for ^^CO(l-O), 
CO (1-0), C^^O(l-O). To complement our study, we have 
also used some data from Ma, Gao, & Wu (2008 in prep.) 
including CS(2-1), HCO+(1-0), HNC(l-O) and HCN(l-O) 
mapping data. These ten sources are marked with • in Ta- 
ble 1 and they were observed with the SEQUOIA on the for- 
mer Five Colleges Radio Astronomy Observatory (FCRAO) 
14m telescope with a 25" on-the-fiy mapping step. The de- 
tailed observation log is listed in Table 2. 

The spectral data were reduced and analyzed with 
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Table 1. Source positions and rms noise levels in spectra (with channel width 50 KHZ). 





RA. 


Dccl. 


Distance 




ccntroid 


CS 2-1 


HCO+ 1-0 


(.18 


O 1-0 


Map Size 


Map Size 


source 


(2000.0) 


(2000.0) 


(kpc) 


Hll? 


(') 


K 


K 




K 


CS 2-l(') 


HCO+ l-O(') 


G123.07-6.31 


00:36:47.51 


-1-63:29:02.08 


2.2 


Hll 


(0,0) 


0.038 


0.068 


0. 


,053 


3x3 


3x3 


•W3(OH) 


02:27:04.69 


-1-61:52:25.54 


1.95 


UCHll 


(0,0) 


0.083 


0.060 


0. 


,038 


7x7 


7x7 


S231 


05:39:12.91 


-1-35:45:54.05 


2.3 




(0,0) 


0.061 


0.084 


0. 


072 


2x2 




S235 


05:40:53.32 


-1-35:41:48.76 


1.6 


Hll 


(0,0) 


0.058 


0.087 


0. 


093 


2x2 




S252A 


06:08:35.41 


-1-20:39:02.96 


1.5 


Hll 


(0,0) 


0.038 


0.081 


0. 


037 


2x2 




• S255 


06:12:53.72 


-1-17:59:22.03 


1.3 


UCHll 


(0,0) 




0.115 


0. 


,066 




7x7 


G19. 61-023 


18:27:38.84 


-11:56:27.44 


4.0 


CHll 


(0,0) 


0.039 


0.034 


0. 


,072 






G20.08-0.13 


18:28:09.88 


-11:28:48.17 


3.4 


UCHll 


(0,0) 


0.040 


0.046 


0. 


,086 






G24.49-0.04 


18:35:05.34 


-07:31:26.98 


3.5 




(0,0) 






0. 


106 






•W44 


18:53:18.50 


-1-01:14:56.65 


3.7 


CHll 


(0,0) 


0.092 


0.104 


0. 


,075 


7x7 


7x7 


S76E 


18:56:10.43 


-1-07:53:14.11 


2.1 


Hll 


(0,0) 


0.040 


0.061 


0. 


,092 


3.5x3 




G35.58-0.03 


18:56:22.57 


-1-02:20:27.75 


3.5 


UCHll 


(0,0) 


0.038 


0.043 


0. 


,076 


2x2 




G35.20-0.74 


18:59:12.73 


-1-01:40:40.75 


3.3 


Hll 


(0,0) 


0.027 


0.027 


0. 


091 






OH43.80-0.13 


19:11:54.27 


-1-09:35:55.38 


2.7 


UCHll 


(0,0) 


0.039 


0.052 


0. 


,078 


2x2 




• S87 


19:46:20.45 


-1-24:35:34.41 


1.9 


UCHll 


(0,0) 


0.111 


0.088 


0. 


,076 


7x7 


7x7 


• S106 


20:27:25.74 


-1-37:22:51.82 


0.6 


UCHll 


(0,0) 


0.070 


0.079 


0. 


,030 


7x7 


7x7 


•W75N 


20:38:36.93 


-1-42:37:37.52 


3.0 


UCHll 


(0,0) 


0.091 


0.071 


0. 


047 


7x7 


7x7 


DR21S 


20:39:00.80 


-1-42:19:29.84 


3.0 


UCHll 


(0,0) 


0.038 


0.076 


0. 


137 


2x2 


2x2 


•W75(OH) 


20:39:01.01 


-1-42:22:49.84 


3.0 




(0,0) 


0.118 


0.090 


0. 


,085 


7x7 


7x7 


• BFSll-B 


21:43:06.68 


-1-66:07:04.11 


2.0 




(0,0) 


0.085 


0.107 


0. 


,065 


7x7 


7x7 


• Ccp-A 


22:56:18.14 


-1-62:01:46.35 


0.7 


UCHll 


(0,0) 


0.117 


0.089 


0. 


,024 


7x7 


7x7 


• NGC7538 


23:13:44.86 


-1-61:26:50.71 


2.8 


UCHll 


(0,0) 


0.232 


0.073 


0. 


,032 


7x7 


7x7 


AFGL4029 


03:02:00.00 


-1-60:33:49.00 


2.2 


UCHll 


(-3,-4) 


0.054 


0.061 


0. 


,063 


2.5x2.5 


2.5x3 


AFGL5142 


05:30:48.00 


-1-33:47:53.80 


1.8 


Hll 


(0,0) 


0.041 


0.073 


0. 


,040 


2x2 


2x2 


S235N 


05:41:00.00 


-1-35:48:04.00 


1.6 




(6,2) 


0.042 


0.068 


0. 


,053 


2x2 


2x2 


G192 


05:58:17.07 


-1-16:31:58.00 


2.0 


UCHll 


(-1,0) 


0.025 


0.052 


0. 


,070 


2x2 




IRAS19410 


19:43:11.00 


-1-23:44:06.00 


2.2 


UCHll 


(0,0) 


0.031 


0.040 


0. 


,028 


3x3 


2x2 


S88 


19:46:43.00 


-1-25:12:14.00 


2.1 


UCHll 


(1,1) 


0.047 


0.072 


0. 


,059 


2x2 


2x2 


IRAS20126 


20:14:25.80 


-1-41:13:33.00 


1.7 




(0,0) 


0.039 


0.114 


0. 


,094 


3.5x3 





NOTE. -positions and distances of sources above the blank line are from Shirley ct al. 2003, except for W3(OH) and S106. The distance of W3{OH) is 
from Xu et al. 2006. The distance of S106 is from Eiroa Elasser & Lahulla 1979; Ridge et al. 2003. For the last eight sources, positions are from Spitzer 
observing center, distances are from Carral et al. 1999, except for IRAS19410. The distance of IRAS19410 is from Xu et al. 2007. Note that S235 has a 
relatively extended distribution, the north region was indi cated by S235N. Ccntroid is the l^CO(l-O) peak position of each source. Columns 7,8,9 are 
rms noise levels twords "'"^CO(l-O) peak position of each source in CS{2-1), HCO"'~(1-0), C"'"^O{l-0) respectively. Symbol ... in the columns 10 and 
11 denoted only single point was observed. 



Table 2. Observed lines and telescope parameters. 



Line 


Data 


Frequency 


Telescope 


FWHM 




a 


Av 






(UT) 


(GHZ) 






(arcsec) 






(kms^l) 


K 


12CO(1-0) 


2005 Nov 


115.271204 


PMO 


14m 


55 


0, 


72 


0.37 






2006 Nov 


115.271204 


PMO 


14m 


55 


0, 


,57 


0.37 




i^CO(l-O) 


2005 Nov 


110.201353 


PMO 


14m 


55 


0, 


,72 


0.11 




2006 Nov 


110.201353 


PMO 


14m 


55 


0, 


,57 


0.11 


200-350 


i*'CO(l-0) 


2005 Nov 


109.782173 


PMO 


14m 


55 


0, 


,72 


0.12 






2006 Nov 


109.782173 


PMO 


14m 


55 


0, 


,57 


0.12 




HCO+{1-0) 


2007 Jun 


89.1885260 


PMO 


14m 


55 


0, 


,62 


0.14 


280-400 




2005 May 


89.1885180 


FCRAO 


14m 


57 


0, 


,53 


0.16 


160-210 


CS(2-1) 


2007 May 


97.9809530 


PMO 


14m 


55 


0, 


.62 


0.13 


160-250 




2005 May 


97.9810070 


FCRAO 


14m 


54 


0, 


,49 


0.15 


120-210 


HNC(l-O) 


2005 May 


90.6635720 


FCRAO 


14m 


57 


0, 


,53 


0.16 


160-210 


HCN(l-O) 


2005 May 


88.6318473 


FCRAO 


14m 


57 


0, 


,53 


0.17 


120-200 



rjmb is the main beam efficiency. 



CLASS and IDL software. Linear baselines were removed 
from all spectra. The rms levels for each spectrum at the 
CO emission peak are also listed in Table 1. The complete 
description of our multiple molecular line observations and 
all data presentation will be given in a forthcoming paper 
(Sun & Gao 2008 in prep.) 

There is little overlap in the samples of previous obser- 
vational search for infall signatures, particularly the map- 
ping studies. Two major surveys of inflow motion observa- 
tions that have some source overlaps with our sample are Wu 



& Evans (2003) and Klaassen & Wilson (2007). There are fif- 
teen sources that overlap between our HCO"^(1-0) mapping 
observations and that of single-pointing HCN(3-2) observa- 
tion of Wu & Evans (2003). Yet, only five sources in our 
sample overlap with the single-pointing HCO^(4-3) obser- 
vation of Klaassen & Wilson (2007) . 
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3 ANALYSIS & RESULTS 

All targets are detected in HCO+(1-0), CS(2-1), C^^O(l- 
0), ^^CO(l-O), ^^CO(l-O), except for G24.49-0.04 owing to 
a wrong frequency setup (a wrong Vlsr was taten). Most 
sources are pretty strong with TX of HCO+(1-0), CS(2-1) 
mainly in a range of 1-4.5 K. But several sources (G19.61- 
023, G20.08-0.13, G35.58-0.03, G35.20-0.74) show fairly 
weak emission (TJ < 0.3 K) in HCO"'^(1-0), and one source 
(G35. 20-0.74) shows quite weak emission (TX < 0.2 K) in 
CS(2-1). Among the 29 dense cores observed, S255 shows a 
very marginal detection in CS omission and S76E has two 
velocity components in both 0^*^0(1-0) and "CO(l-O) and 
is difHcult to determine the systematic velocity. Therefore, 
these sources (G24.49-0.04, S255, S76E) are excluded from 
the following analysis. We only give the results derived from 
CS(2-1) and HCO^(l-O) as optically thick spectral lines to- 
wards 26 sources below since only ten sources have both 
HCN(l-O) and HNC(l-O) observations. 



3.1 Blue profile identification 

To quantify the blue asymmetry of a line, we use an asym- 
metry parameter 5Y (pioneered by Mardones et al. 1997) de- 
fined as the difference between the peak velocities of an opti- 
cally thick line V(thick) and an optically thin line V(thin) in 
units of the optically thin line FWHM (Full Width at Half 
Maximum) dV(thin): SV = ^ Mardones et 

al. (1997) adopted a criterion 5V < —0.25 to indicate blue 
asymmetry and SV > 0.25 to indicate red asymmetry. They 
chose the threshold value 0.25 at about 5 times the typical 
rms error in 5V in order to exclude the random contributions 
(Mardones et al. 1997). We use the same criterion for con- 
sistency with their analysis. Because of the asymmetric pro- 
files for many of optical-thick lines, they cannot be fitted by 
Gaussian function, the peaks of the line profiles of CS(2-1), 
HCO+(1-0), and HNC(l-O) are identified with the highest 
flux by cursor. The uncertainties of the velocities and the an- 
tenna temperatures at these peaks were estimated by several 
measurements and the baseline fitting noises. The optically 
thin line (C^^O) shows an ideal Gaussian profile and is fit- 
ted by the Gaussian function to obtain the peak velocity 
and FWHM. All these parameters towards the "CO(l-O) 
peak position are listed in Table 3 which includes the peak 
velocities of CS(2-1), HCO+(l - 0), C^®O(l-0), the FWHM 
of C^*O(l-0), the calculated asymmetry parameter SV of 
CS(2-1) and HCO+(1-0), and the uncertainties of these pa- 
rameters (in the parentheses) for all sources. 

The blue profile caused by infall motion also requires 
T^iB) /T'^[R) > 1 -f (7 (i.e., the blue peak is stronger than 
the red peak, hereafter Tb/Tr = Tl{B)/Tl{R)). These ra- 
tios of CS(2-1) and HCO+(1-0) towards the ^^CO(l-O) peak 
position are available in Table 3. The peak antenna temper- 
atures of CS(2-1) and HCO"'"(1-0) towards each source arc 
also listed in table 3 columns 8 and 10, respectively. Based 
on the methods that wc used to derive the ratio of Tb/Tr, 
the sources can be classified into three groups: (1) sources 
show obvious two peaks, we can accurately measure the ra- 
tio of Tb/Tr; (2) sources with asymmetry profiles but do not 
show obviously two peaks (both blue and red) , thus only one 
peak (either blue or red velocity offset) can be identified. We 
simply chose the same velocity offset to the systematic veloc- 
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Figure 1. Distribution of <5V (left, the line asymmetry parame- 
ter) and log[TX(B)/T;j^(R)] (right, ratio of blue to red peak inten- 
sity). The top and bottom panels show results from CS(2-1) and 
HCO+(1-0) respectively. 



ity as the other peak, so Tb/Tr can not be exactly derived. 
We mark "~" before such sources; (3) sources almost don't 

show any asymmetry, neither the rod nor blue component 
can be identified, so the intensity ratio Tr/Tr was flagged 
with 

In our sample, four blue profiles were identified in 
CS(2-1), fourteen blue profiles were identified by HC0"'"(1- 

0), three sources show blue profiles both in CS(2-1) and 
HCO^(l-O) spectra. Within the error bars, the parame- 
ters (<5V, Tb/Tr) of several sources (W44, G192, G35.20- 
0.74,S88) are approaching the criterion, so they are possible 
candidates of the blue profiles or red profiles and are marked 
with "?" in Table 3. Note that S87 and G35.20-0.74 show 
different asymmetry between CS(2-1) and HCO+(1-0). S87 
also shows different asymmetry at different positions. At off- 
sets (0", 25"), (0", -25"), (25", 0"), (-25", 0"), it shows red 
asymmetry, blue asynnnetry, no asymmetry, blue asymme- 
try in HCO^(l-O) and red asymmetry, no asymmetry, red 
asymmetry, red cisymmetry in CS(2-1), respectively. These 
might suggest a complicated kinematics in the mapped re- 
gions in S87 (c.f. Xue & Wu 2008, where they claimed that 
the asymmetric line profiles in S87 should result from two 
clouds at slightly different velocities). Figure 1 presents the 
statistics of the asymmetry parameters measured towards 
the '^^CO(l-O) peak position of the 26 sources observed in 
both HCO+(1-0) and CS(2-1) mentioned before. It's clear 
that the blue profile predominance was shown, particularly 
in HCO'''(1-0). Because infall is the only process that would 
produce consistently the blue profile, its predominance in- 
dicates that the infall motions are really present in mas- 
sive star-forming regions. A comparison of the asymmetry 
parameter S\f in CS(2-1) and HCO^(l-O) is shown in Fig- 
ure 2, which shows weak correlation between 5v in the two 
molecular lines. As noted above, G35.20-0.74 and S87 show 
different asymmetry between CS(2-1) and HCO'''(1-0) and 
also greatly affect the correlation between (5V of the two 
lines. Generally, there are more sources with blue asymmet- 
ric profiles in HCO+(1-0) Une than those in CS(2-1). We will 
further discuss the difference between the two lines in S4. 
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Table 3. The Derived line parameters and profiles of the observed sources. Quantities in parentheses give the uncertainties in units of 0.01. 



source ^ 




Vthick 




AV 


Sv 


i5v 




Tb/Tr 




Tb/Tr 


profile ^ 




CS 2-1 


HCO+ 1-0 


C^^O 1-0 


C^^O 1-0 


CS 2-1 


HCO+ 1-0 


CS 2-1 


CS 2-1 


HCO+ 1-0 


HCO+ 1-0 






km s ^ 


km 


km s ^ 


km 






K 




K 






*gl23.07-6.31 


-18.04(04) 


-18.35(01) 


-17.55(02) 


1.80(04) 


-0.27(03) 


-0.44(02) 


1.52 


1.54 


2.95 


1.79 


B,B,B 


*W3(OH) 


-47.53(06) 


-48.76(07) 


-47.15(02) 


4.16(05) 


-0.09(02) 


-0.39(02) 


3.73 


~1.08 


2.87 


2.34 


N,B,B 


S231 


-17.14(03) 


-17.33(01) 


-16.76(03) 


3.66(07) 


-0.10(01) 


-0.16(01) 


1.31 


~1.08 


2.18 


~1.09 


N,N,R 


S235 


-16.83(09) 


-17.30(04) 


-16.85(04) 


2.38(09) 


0.01(04) 


-0.19(02) 


2.12 




2.52 


~1.47 


N,N,N 


S252A 


8.95(05) 


8.60(10) 


8.53(01) 


2.44(03) 


0.17(02) 


0.03(04) 


1.76 


~0.91 


1.56 




N,N,R 


G19.61-0.23 


42.79(12) 


40.31(06) 


42.75(12) 


6.61(34) 


0.01(03) 


-0.37(03) 


0.62 




0.32 


~1.58 


N,B,N 


G20.08-0.13 


41.79(14) 


41.29(07) 


42.13(13) 


4.83(39) 


-0.07(04) 


-0.17(03) 


0.53 




0.29 


~1.53 


N,N,~ 


oW44 


58.02(13) 


56.51(07) 


57.76(02) 


5.20(04) 


0.05(03) 


-0.24(01) 


3.27 




2.96 


~4.96 


N,B?B 


G35.58-0.03 


53.29(04) 


54.00(15) 


53.01(02) 


6.16(06) 


0.05(01) 


0.13(01) 


1.01 


~0.94 


0.34 


~0.39 


N,N,N 


G35.20-0.74 


37.50(04) 


36.25(09) 


37.09(05) 


1.35(09) 


0.30(05) 


-0.62(09) 


0.16 


~0.79 


0.17 


~1.05 


R?B?~ 


OH43.80-0.13 


43.88(06) 


43.82(12) 


44.41(12) 


5.30(21) 


-0.10(03) 


-0.11(03) 


0.72 


~1.19 


0.80 


~1.20 


N,N,~ 


S87 


23.35(06) 


21.20(07) 


22.42(09) 


3.39(16) 


0.27(03) 


-0.36(04) 


2.15 


0.47 


3.28 


1.26 


R,B,~ 


0SIO6 


-1.47(12) 


-1.85(07) 


-1.01(01) 


2.47(03) 


-0.19(05) 


-0.34(03) 


1.74 


~1.58 


1.63 


~1.33 


N,B,~ 


W75N 


9.25(07) 


9.55(05) 


9.69(02) 


3.47(05) 


-0.13(02) 


-0.04(02) 


3.30 


~1.00 


2.97 




N,N,R 


DR21S 


-2.59(03) 


-4.47(05) 


-2.78(04) 


3.12(09) 


0.06(02) 


-0.54(03) 


2.88 




2.41 


1.35 


N,B,B 


*W75(OH) 


-4.28(07) 


-5.03(07) 


-3.06(03) 


2.97(07) 


-0.41(03) 


-0.66(03) 


4.36 


1.88 


4.11 


2.00 


B,B,B 


BFSll-B 


-10.05(05) 


-10.25(05) 


-10.22(02) 


1.58(08) 


0.11(03) 


-0.02(03) 


1.50 


~1.25 


2.35 




N,N,~ 


*Ccp-A 


-11.44(16) 


-12.71(10) 


-10.62(01) 


3.56(03) 


-0.23(05) 


-0.59(03) 


2.43 


1.19 


2.29 


2.29 


N,B,B 


*NGC7538 


-57.36(14) 


-57.81(21) 


-56.00(03) 


5.40(07) 


-0.25(03) 


-0.34(04) 


5.16 


~1.48 


3.96 


~3.42 


B,B,~ 


AFGL4029 


-38.35(01) 


-38.49(03) 


-38.08(04) 


2.19(10) 


-0.12(02) 


-0.19(02) 


1.89 


~1.19 


1.82 


~1.22 


N,N,~ 


AFGL5142 


-3.51(02) 


-3.84(02) 


-3.37(02) 


2.75(04) 


-0.05(01) 


-0.17(01) 


3.67 




3.71 


~1.31 


N,N,~ 


*S235N 


-20.23(10) 


-20.83(04) 


-19.62(05) 


2.83(11) 


-0.22(04) 


-0.43(03) 


1.08 


~1.28 


1.98 


2.68 


N,B,~ 


G192 


5.61(06) 


6.28(05) 


6.18(05) 


2.10(10) 


-0.27(04) 


0.05(03) 


0.31 


~1.63 


0.70 


1.09 


B?N,~ 


IRAS19410 


22.54(01) 


21.50(04) 


22.57(01) 


2.35(02) 


-0.01(01) 


-0.46(02) 


1.82 




1.50 


1.06 


N,B,B 


S88 


21.90(06) 


21.38(11) 


22.04(03) 


2.69(06) 


-0.05(02) 


-0.25(04) 


0.91 




2.43 


2.18 


N,B?R 


IRAS20126 


-3.88(03) 


-3.65(05) 


-3.58(06) 


2.45(16) 


-0.12(03) 


-0.03(03) 


1.06 


~1.16 


2.24 




N,N,~ 



" Sources marked with o and * are considered to be probably and strong infall candidates, respectively. 

^ Profile is judged from our CS(2-1), HCO+(1-0) observations and HCN(3-2) observation of Wu & Evans (2003), respectively. B denotes 
blue profile, R denotes red profile, N denotes neither blue nor red, ~ denotes no HCN(3-2) data 




-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 

^ ^ HCO (1-0) 



Figure 2. Comparison of the measured asymmetry in HC0"'"(1- 
0) and CS(2-1). The dashed lines mark \5V\ = 0.25 and the line 
of equality is shown in solid. 

3.2 Infall candidates 

The mapping observations of multiple transitions could offer 
detailed investigation to provide strong evidence on whether 
the core is an infall candidate. It has been widely accepted 
that the infall candidates must show obvious blue profile at 
least in one optically thick line, no red profile in other infall 
tracers and no spacial difference in the mapping observa- 
tion. Based on all of these, we have found three probable 
infaU candidates (S106, W44, IRAS19410) and six strong 



infaU candidates (G123.07-6.31, W75(OH), S235N, CEP- 
A, W3(0H), NGC7538). The strong infaU candidates are 
marked in Table 3 and Figure 2. 

Figures 3-8 show the HCO^(l-O) mapping observations 
towards the six most strong infall candidates (Figs. 4-7 are 
presented with only central parts of the HCO"'' (1-0) maps). 
Most of the data show obvious asymmetric profiles and all of 
spectral lines observed towards the CO (1-0) emission peak 
position. The six infall candidates show consistent blue pro- 
file in at least one optically thick spectral line. Most of them 
show high velocity wing emissions both in HCO"'' (1-0) and 
^■^ CO (1-0). These might suggest that besides infall motions 
outflows are also present in the infall candidates. Klaassen 
& Wilson (2007) also found that seven out of their nine in- 
fall candidates show SiO detections. All these confirm that 
infiow was often accompanied by outflow in the process of 
massive star forming process. In these strong infall candi- 
dates a large scale blue asymmetry was present. From our 
maps we estimated the extent of the infall signature to be 
up to 2' X 2' (at least 0.66 pc for G123. 07-6.31, 0.59 pc for 
W3(OH), 1.03 pc for W75(OH), 0.21 pc for CEP-A, 0.84 pc 
for NGC7538, and 0.34 pc for S235N). 

We here choose G123.07-6.31 as a typical example of 
infall candidate to illustrate our observations (Figure 3). In 
Figure 3a, a signature of extended inward motions were de- 
tected in all the mapping region of our HCO^(l-O) obser- 
vation, detailed analysis was needed to better understand 
the spatial variation of the infall speed of the core. Figure 
3b shows all of molecular lines observed towards the maser 
site (0,0), and outer region of the core (l',-l')- Obvious blue 
asymmetry is presented in both CS(2-1) and HCO''' (1-0). 
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Along Adecl. =0 the position-velocity (P-V) diagram of 
HCO^(l-O) shows strong blue asymmetry (Figure 3c top) 
while the P-V diagram of C^*O(l-0) doesn't show any asym- 
metry (Figure 3c bottom). 

A noticeable feature in Figure 9 is that the HCO"''(1-0) 
lines generally have higher intensity ratios Tb/Tr than that 
in CS(2-1) for the six infall candidates; typically~2.2 for 
HCO+(1-0) and ~1.4 for CS(2-1). The difference is signifi- 
cantly larger than any possible errors produced in determin- 
ing Tb/Tr. a simple analytic model for collapse (Myers et 
al. 1996) predicted that the ratio of blue to red component 
(Tb/Tr) increases with infall speed (Vi„), which assumes 
that the cloud consists of two uniform paralleled compo- 
nents without rotation and the infall speed of the gas is less 
than the velocity dispersion a. Thus, it is likely that HCO-I- 
and CS trace different infall speeds given their different in- 
tensity ratios. This also suggests that the infall speeds de- 
rived from HCO+(1-0) and CS(2-1) might be different. Of 
course, we do not expect the derived infall speeds to be same 
because different infall tracers might trace different spatial 
components of the clouds along the line of sight. From the 
mapping observations, it is obvious that the infall speed is 
also different in different spatial components perpendicular 
to the line of sight. Thus, these two lines do not necessarily 
provide duplicate information. Rotation might cause blue 
asymmetric line profile at one side of the rotation axis, but 
at the same time the red asymmetry could appear at the 
other side of the rotation axis. Based on this argument, we 
could exclude the rotation as the main contributor. We did 
not do any model fit towards these sources in this paper, 
some modelling efforts will need to be developed in order to 
derive the infall motion and mass infall rate in future work. 



4 DISCUSSION 

As presented in the previous section, an obvious asymmetry 
parameter difference was shown between HCO^(l-O) and 
CS(2-1) in Figs. 2 and 9. HCO+(1-0) line seems to be a 
better tracer for the study of infall motion compared with 
CS(2-1). In general, the characteristic blue profile will ap- 
pear only if the molecular line tracer has a suitable optical 
depth and critical density. Both HCO+(1-0) and CS(2-1) are 
optically thick lines, but HCO'''(1-0) is more opaque than 
CS(2-1) (Evans 2003). The critical density of the two lines 
is also different. At Tfe = 100K, the critical density of CS(2-1) 
is 3.9 X 10^ cm~^ (Evans 1999), while the critical density of 
HCO+(1-0) is 1.9 X 10^ cm"^ (Evans 1999; cf. 5.0 x 10'' cm"^ 
Fuller et al. 2005). A natural explanation for the asymme- 
try parameter difference between CS(2-1) and HCO'''(1-0) is 
that infall motions are more easily detected in the lower den- 
sity outer region of the dense core. In the denser regions, the 
infall motions might be significantly decelerated by thermal 
pressures, outflows, and other feedbacks from the massive 
young objects. 

Nevertheless, HCN(3-2) has a much higher critical den- 
sity (6.8 X 10^ cm"^ at Tfe = 100 K) than that of CS(2-1), 
yet it still shows a high inflow motion detection rate (Wu & 
Evans 2003). In the flfteen sources that overlap between our 
survey and that of Wu & Evans (2003), the profiles derived 
from 1100"*" (1-0) are more consistent with that derived from 
HCN(3-2) than that of CS(2-1) (see Table 3.). The difference 
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Figure 3. G123.07-6.31 (a)HCO+(l-0) map grid. (b)Spectral to- 
wards (0,0) (left), and (l',-l') (right) the dashed line marked 
the source systemic velocity which was derived from gauss fit 
of Ci*O(l-0). (c)Position-velocity dia grams of optically thick 
HCO+(1-0) line (top) and optically thin 0^*0(1-0) line (bot- 
tom), contour levels are 0.4K (> Scr), 0.7K, l.OK, 1.3K, 1.7K, 
2. IK, 2.5K, 2.9K and 0.3K (> Scr), 0.5K, 0.7K, 0.8K, 0.9K, re- 
spectively. 
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Figure 4. W3{OH), the top and bottom panels present the center 
parts that show obvious asymmetric profiles of HCO+(1-0) map 
grid and all of spectral lines observed towards (0,0), respectively. 
In the top panel, the velocity scale ranges from -62 km to 
-32 km s-^ the same as the bottom panel, the temperature scale 
ranges from -0.5 K to 4.4 K. 



between the optical depth of the two molecular lines might 
be the more important reason than the diflterence in the 
critical density of the two lines. A recent model about blue 
asymmetry in low-mass star-forming region suggests that 
the line asymmetry is strongly dependent on the abundance 
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Figure 5. W75{OH), the top and bottom panels present the cen- 
ter parts that show obvious asymmetric profiles of HCO+(1-0) 
map grid and all of spectral lines observed towards (0,0), respec- 
tively. In the top panel, the velocity scale ranges from -20 km 
to 15 km s^^ the same as the bottom panel, the temperature 
scale ranges from -0.5 K to 4.4 K. 



distribution (Tsamis et al. 2008), which is proportional to 
the optical depth. 

In the fifteen overlap sources, more blue and red profiles 
were presented by Wu & Evans (2003) than that of ours. 
Their HCN(3-2) observations were carried out in a single 
pointing toward dense cores using the CSO 10.4 m telescope 
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Figure 6. CEP-A, the top and bottom panels present the center 
parts that show obvious asymmetric profiles of HCO+{1-0) map 
grid and all of spectral lines observed towards (0,0), respectively. 
In the top panel, the velocity scale ranges from -40 km s~^ to 
20 km s-l the same as the bottom panel, the temperature scale 
ranges from -0.5 K to 4.4 K. 
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Figure 7. NGC7538, the top and bottom panels present the cen- 
ter parts that show obvious asymmetric profiles of HCO+(1-0) 
map grid and all of spectral lines observed towards (0,0), respec- 
tively. In the top panel, the velocity scale ranges from -80 km 
to -32 km s~^ the same as the bottom panel, the temperature 
scale ranges from -0.5 K to 4.4 K. 



with resolution about 28" . The higher spatial resolution en- 
ables them to better resolve the inner small regions of in- 
ward motions. In particular, G19.61-0.23 shows blue profile 
in both HCO+(1-0) of our results and HCO+(4-3) of the 
JCMT observations with resolution about 15" (Klaassen & 
Wilson. 2007), but no obvious asymmetry in HCN(3-2) (Wu 
& Evans 2003). S88 shows blue profile in HCO+(1-0) over a 



2' X 2' mapping region in our observations whereas it shows 
red profile in a single pointing HCN(3-2) observation of Wu 
& Evans (2003), yet it shows no obvious asymmetry in sin- 
gle pointing HCO'''(4-3) observation of Klaassen & Wilson 
(2007). In a recent survey, Fuller et al. (2005) also found 
that the detection rate of infall motion decreases with in- 
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Figure 8. S235N, the top and bottom panels present HCO+(1-0) 
map grid and all of spectral lines observed towards (6', 2'), respec- 
tively. In the top panel, the velocity scale ranges from -40 km 
to km s~^ the same as the bottom panel, the temperature scale 
ranges from -0.5 K to 2 K. 



creasing source distance. In short, the source distance, half 
power beam width (HPBW) of the molecular transition line 
probes, and the size scale of infall region could all affect the 
infall detection rate. 

Even though our observations (resolution ~ 1') did not 
have the resolving power to probe a linear scale <0.16 pc, 
the mapping observation of multiple optically thick molec- 
ular line could still be useful to give evidence to large 
scale infall signatures in dense cores. Six infall candidates 
were identified in our mapping observations, in which large 



CS(2-1) ■ 
HCO*(1-0)^ 



I 



Figure 9. Intensity ratios Tb/Tr in HCO+(1-0) and CS(2-1) for 
the six infall candidates. 



scale (>0.2pc) infall signatures were detected. Among them, 
NGC7538 has been previously well studied by Sandell et al. 
(2005). Their high resolution BIMA HCO+(1-0) observation 
suggests that an embedded protostar embedded in the core 
is still in a phase of active accretion. All these might sug- 
gest that large scale infall motions are closely related to the 
kinematics of the inner materials of the cores. A large scale 
infall might play a certain role in providing some materials 
for the innermost accretion. 

Actually, outflow makes an important contribution to 
the line wing emissions. But high quality spectra of much 
higher signal-to-noise ratios are needed to make better 
judgement in line wings. Following previous similar studies 
and our observational results, we here believe that the line 
asymmetry is mainly contributed by infall motion. And to a 
certain extent, large scale infall signature at least provides 
an indirect evidence for the inner accretion. Protostellar ac- 
cretion could be the main physical process responsible for 
massive star formation. Some models suggest that the accre- 
tion must halt before a visible UC Hll phase (eg. Garay & 
Lizano 1999), while other models suggest that an ionized ac- 
cretion can continue through an Hll region (eg. Keto 2003). 
Four out of our six infall candidates have already formed 
UC Hii region (W3(OH), CEP-A, NGC7538) and Hii region 
(G123.07-6.31), yet they still show large scale infall signa- 
tures. Although we don't know whether infall motions halt 
at smaller scale of the cores, these examples at least at large 
scale appeared to be consistent with the ionized accretion 
theory of Keto (2003). 

Comparing with other surveys, our observations show 
only two (8%) red asymmetry detection in CS(2-1) and no 
red asymmetry detection in HCO+(1-0) (vs. 21% in HCN(3- 
2) observation of Wu & Evans (2003); 16% in HCO+(l- 
0) observation of Fuller et al. (2005); 10% in HCO+(1-0) 
of Wu et al. (2008)). The likely reason is beam dilution 
and/or the fact that the samples of all surveys until now 
are still rather limited (typically ~30, e.g., 28 sources of Wu 
& Evans (2003); 23 sources of Klaassen & Wilson (2007), 
vs. 77 sources of Fuller et al. (2005)). Of course, these are 
just some face values and the exact percentage should be 
affected by resolution. But it is rather arbitrary to scale the 
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resolution in comparing to other observations. S87 is such 
an interesting source which shows different asymmetry in 
CS(2-1) and HCO+(1-0) at different offsets (so might be 
G35.20-0.74). With our telescope resolving power, it's hard 
to give a reasonable explanation for the different asymmetry 
in different positions. The possible explanation also include 
cloud-cloud coUision (e.g. S87, Xue & Wu 2008). Oscillat- 
ing was also suggested in some low mass cores, e.g., B68 
(Redman et al. 2006) and perhaps this might be possible 
in these massive cores as well, ffigli resolution observations 
are required towards these sources in order to learn more 
details about the kinematics of massive star-forming cores. 
Unlike blue profile, the origin of red profile is still unclear. 
A simple explanation is that tfie red profile could be caused 
by outward motions. 



5 SUMMARY 

We performed a systematic multiple transition 3 mm molec- 
ular line single-pointing and mapping survey towards 29 
massive star-forming cores to search for the signatures of 
the inward motions. Up to seven different transitions, op- 
tically thick lines HCO+(1-0), CS(2-1), HNC(l-O), HCN(1- 
0), ^^00(1-0) and optically thin lines C^^O(l-O), ^•'CO(l-O) 
were observed towards each source. The major results found 
from these observations include the following: 

(i) The normalized velocity differences {SVcs, SV^cq+) 
between the peak velocities of optically thick lines and opti- 
cally thin C^®O(l-0) for each source were derived. We found 
a significant difference in the incidence of blue-shifted line 
asymmetry between CS(2-1) and HCO+(1-0). The HCO+(l- 
0) line shows the highest occurrence of obvious asymmetric 
feature. The optical depth difference between the two molec- 
ular lines may be responsible for the anomaly. HCO^(l-O) 
line appears to be one of the best inward motion tracers in 
massive cores whereas CS(2-1) is not a sensitive tracer to 
signify the infall motion. 

(ii) Blue profile dominance appeared in both CS(2-1) and 
HCO^(l-O) in most of our sample. This implies the predom- 
inance of the inward motion in the massive star-forming 
cores, which is similar to that in low mass star forming 
regions. All those may suggest that protostellar accretion 
could be the main physical process responsible for massive 
star formation. 

(iii) The mapping observations of multiple line transitions 
enable us to identify six strong infall candidates (G123.07- 
6.31, W75(OH), S235N, CEP-A, W3(OH), NGC7538). Sig- 
natures of extended inward motions are presented in these 
massive star-forming cores with possibly strong infall. The 
infall signature is extended up to a linear scale > 0.2 pc. 
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APPENDIX A: IRAC AND CO LINE 
COMPARISON IN S235 

The recently available high-resolution and very sensitive 
mid-infrared and far-infrared images of the Spitzer telescope 
can provide a new insight in directly locating the distribu- 
tion of dust obscured massive stars in the dense cores. As an 
example for some well studied dense core sources, we here 
compare our molecular line observations with the Spitzer 
IRAC images in S235. The Spitzer IRAC data were retrieved 
from Spitzer Science Centeo IRAC 4 bands data of S235 
(PI G. Fazio, PID=201) covering~0.25 deg^ reveal a com- 
plex star-forming region. 

The S235 star-forming region is a well known Hll region 
that has been widely studied. It has been identified as a 
water maser (Lo, Burke, & Haschick 1975) with infrared 
(e.g. Krassner, Pipher & Sharpless 1979; FeUi et al. 1997; 
Allen et al. 2005), radio continuum (Felh et al. 1997, 2006), 
dust continuum emission (Lee, Young & Shirley 2002), and 
many molecular line emission (e.g. Evans et al. 1981; Felli 
et al. 1997; Lee, Young, & Shirley 2002). 

Our CO maps show that three CO (1-0) cores are well 
matched with three infrared emission regions in Figure Al. 
We observed several points in these three cores in CS(2-1) 
and HCO+(1-0), only HCO+(1-0) show obvious blue pro- 
file around position (6', 2') marked dash open circle. Class I 
and class II sources in S235 were identified by Allen et al. 
2005, and their results show that few YSOs are distributed 
in the North- West extended core and only two class I sources 
are distributed around position (6', 2'), whereas many more 
YSOs are distributed in the other two cores. All these may 
imply that the region around position (6', 2') is still on the 
very early stage of star formation. 
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Figure Al. The l3cO(l-0) contour map overlaid on the Spitzer 
IRAC 4.5 fim mosaic of S235. Contour levels start at 30% of the 
peak ^^CO(l-O) integrated intensity and increase by 10% at each 
interval. 



^ |http: / /ssc.spitzer. caltech.edu 



